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We report on a possibility of measuring parity violating effects in radiative neutrino pair emission 
from metastable atoms; asymmetric angular distribution of emitted photons from oriented atoms 
and emergent circular polarization. Their observation, along with the continuous photon energy 
q ■ spectrum which has 6 thresholds, may be interpreted as events being a combined weak and QED 
^ ■ process, emission of jUiUj in the final state. The method may greatly help to perform neutrino mass 

• • ! spectroscopy using atoms, a systematic determination of the neutrino mass matrix. 
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Introduction Despite exciting experimental developments in neutrino physics the neutrino 
is still a very mysterious particle. Do neutrinos obey the same type of Dirac equation as all other 
charged leptons and quarks, or do they obey two-component Majorana equation in which particle is 
indistinguishable from its antiparticle ? How small is the lightest neutrino mass ? These questions are 
linked to the cosmological scenario of leptogenesis [1] which ultimately aims at a deep understanding 
of the matter- antimatter imbalance in our universe. 

Radiative neutrino pair emission (RNPE) [2], [3] \i)a — > \I)a + 7 + VxVj (v>i being the i-th mass 
eigenstate ordered in mass by < m 1 < m 2 < m 3 ) from a metastable atomic state \i)a, a process 
yet to be discovered, gives the opportunity of systematically measuring the neutrino mass matrix 
(described by 9 real parameters) and the mass type, Majorana vs Dirac. In the present work we 
calculate the magnitude of observables that evidently indicates parity violation, providing, along 
with a missing energy suggestive of invisible neutrino emission, a fundamental tool of unambiguous 
detection of the weak RNPE process. The RNPE process furthermore gives critical information on 
three neutrino masses and the unknown mixing angle #13 . This way one may advance a step towards 
verifying important ingredients of leptogenesis theory. 

The first experimental discovery of parity violation was made by measuring emitted electron 
correlation with polarization axis of target in the (3 and fi decay [I]. The measured electron in 
the P and [i decay is replaced here in RNPE by the photon. In RNPE the circular polarization 
of emitted photon is another good quantity for verification of parity violation. Mass thresholds of 
invisible neutrino pair emission of mass eigenstates z/jZ/j are measurable in atomic experiments by a 
spectrometer with good energy resolution, while in j3 and /1 decays only flavor eigenstates u e , are 
observable with poor energy resolution. 

In order to increase feasibility of actual experiments one needs enhancement of RNPE rates, for 
which we use macro-coherence of target atoms [i3]. Measurement of parity violation helps much to 
reduce QED backgrounds even for smaller rates, which thus becomes important together with other 
means of the background rejection [5]. 

Prepared states In discussions that follow it is essential to have the momentum conservation 
k + pi + P2 = among 3 emitted light particles, k the momentum of photon, pi those of neutrinos, 
taking infinitely heavy atom, along with the usual energy conservation. This way the macro-coherent 
decay is much similar, in kinematics, to the ordinary decay of a single relativistic particle of mass A 
where A = E,i — Ef is the energy difference of initial and final atomic states. For the momentum- 
conserving three-body decay of RNPE the energy spectrum of photon has 6 different thresholds which 
occur at photon energy of A/2 — (m 8 + rrij) 2 / (2A). The threshold locations thus determine neutrino 
pair masses rrii + rrij. It turns out that the threshold rise of single photon spectrum is sensitive to 
neutrino mixing angles such as the unknown #13 (see below for more details). How the momentum 
conservation holds for macro-coherent phenomena [3] shall be explained later. 

Atoms we consider belong to a class of A— type three level system, \i)a, \ u )a, \I)a with energy 
relation E n > Ei > Ff. RNPE occurs in second order perturbation, a combination of weak and 
QED process, and has the transition operator of the structure Gpd ■ ES e ■ j u / (A n j — uS) , where 
the electron spin operator S e arises from the neutrino pair emission vertex for atomic transition 
\i)A — * \ n )A, multiplied by the neutrino current j u , a bilinear form of two plane wave neutrinos, 
and d ■ E for \u)a —> |/)a is the usual QED El (or Ml /2 ■ B) vertex. Here A n f = E n — Ef is the 
energy difference between two atomic levels, and 10 is the photon energy. We may classify the weak 
vertex as Gamow- Teller (GT) type, since it involves the electron spin operator. The Fermi type 
operator appears as a small relativistic correction, and it usually gives much smaller rates than GT, 
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GT/Fermi interference being of order 10~ 4 less usually. The transition operator for v^Vj emission 
is multiplied by £y = U^Uej — <$y/2, with U ei the neutrino mixing matrix [6]. For simplicity we 
assume that the Majorana phases are absent, thus U e \ = C12C13 , U e 2 = S12C13 , U e % = S13 , with 
Cij = cos% , Sij = sin 9ij ([7] for an effect of Majorana phases). 

Two energy differences of atomic levels, Aj/ = Ei — Ef = A and A n j = E n — Ei, are involved 
in the transition matrix element of RNPE; the larger A and the smaller A n j, the larger the rate is. 
Low lying excited levels of rare gas atoms are among the best candidate targets. For instance, an 
example of Xe gives a good candidate of states; \i)a = 5p 5 ( 2 P3/2)6s 2 [3/2] 2 of total electron angular 
momentum J = 2, metastable with lifetime O[40] sec, 1/)^ = 5p 61 S the ground state of Xe atom, 
and \tl)a = 5p 5 ( 2 P3/2)6s 2 [3/2]i of J = 1. These low lying excited states are well described by two 
"particle" picture of electron and hole; (6s)(5p) + . The system is nearly degenerate, A ni <C A, giving 
a large RNPE rate. In isolated environment A n j ~ 0.12 eV and A ~ 8.4 eV for Xe. 

Preparation of atomic polarization in excited states is not difficult in optical experiments using 
atoms, made possible by tuning excitation laser frequency, if a magnetic field (its direction taken to 
|| z-axis) is applied to distinguish Zeeman sub levels. Both for Xe atoms of even (having no hyperfine 
split levels) and odd nuclei (both HFS and Zeeman splitting present) we may use the total angular 
momentum and its component along the magnetic field F, M to distinguish these energy levels. 

A complication that arises to obtain large PV effects is that one needs a nearly equal mixture 
of states of magnetic quantum numbers M, different by ±1. This is made evident below in explicit 
computations. There may be a variety of methods to realize this mixture, and we consider here a 
simultaneous application of tilted electric field. As is shown in Figure 1 for odd 131 Xe of nuclear spin 
/ = 3/2 [8], the crossing of different hyperfine levels occurs at ~ 1.5 kG for 1 7/2, —1/2) and 1 5/2, 1/2), 
creating degeneracy (we consulted [9] for this computation). Zeeman splitting due to the magnetic 
field is the first order effect of the field, while the Stark effect due to the electric field is of the second 
order in field. Presence of the weak electric field removes the degeneracy at the level crossing point, 
generating a state of nearly equal mixture of states of the form, cosr/|P, M) + sinr/|P — 1, M + 1) 
of tan f] ~ 7r/4. There are actually two nearly degenerate states, this one and the other orthogonal 
one. The orthogonal one gives PV effects exactly opposite in sign, hence one needs to separate two 
contributions. This is experimentally possible, for instance, by adiabatically changing the magnetic 
field strength and choosing either of the two branches of level repulsion. Parity violating effects then 
appear as oc cos?] sin 77. 

Relevant electronic transitions are of the type J = 2-^J = 1— > J = (note AF = AJ). 
Suppose that the magnetic quantum number along z-axis is denoted by M. We wish to obtain 
mixture of states, M = —1/2 and 1/2 states for finite PV effects. The electron hole in the inner 
core (5p) + of the initial state has J = 3/2 for relevant J = 2 —>■ 1 transition. For initial M — ±1/2 
mixture one may consider a hole in p_ and p z whose wave functions are oc (x — iy)/\/2 and oc z. 
Suppose that a weak Stark field E s is applied to a tilted direction in (yz) plane, with a tilted angle ip 
from z-axis. Second order Stark effects caused by Hamiltonian E s (z cos + ysinip) then leads to an 
off-diagonal mixing of the type, 6s(5p_) + — > (6s) 2 — > 6s(5p z ) + giving a mixture of states, 6s(5p±) + 
and 6s(5p z ) + . The matrix element for this is oc cosip sin ip, thus is maximal at tp — tt/4. In the rest 
of numerical estimates we assume ip = tt/4. 

RNPE with invisible neutrino pair For illustration we consider 131 Xe atomic transitions 
of the type, 
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Figure 1: Zeeman levels of odd 131 Xe atom. 



Summation over 4 final states is inevitable in a simplest experimental setup that does not measure 
final atomic states. 

Neutrinos are difficult to detect, but one may measure energy, emission angle, and circular po- 
larization of the accompanying photon in RNPE, and derive event rates of specified photon energy- 
momentum, k, and helicity /i 7 , in the form of continuous spectrum. On the other hand, we sum 
over non-observable variables; neutrino helicities hi and their momenta Pi of the pair. The neutrino 
helicity summation in rates j2] gives bilnear terms of neutrino momenta; Kf- = 

£ jWm) ] = (-^s M + i)<% + -^riplti + M - ■ n) . (2) 

ftl,ft2 

The case 5m = 1 is applied for Majorana neutrinos, and 5m = for Dirac neutrinos. The sub- 
sequent neutrino momentum integration (with E{ = \J pf + rnf neutrino energy), with J dP v = 
J d 3 pid 3 p 2 5 3 (k + f>i + p 2 )5(A — uj — Ei — E 2 )Kf- , gives a 1-dim. neutrino energy integral f dE\, 
which can be done analytically. 

The result of neutrino helicity and momentum summation gives quadratic forms of photon mo- 
mentum, kikj and 5ijk 2 x functions of photon energy \k\ = uj. These quadratic forms are to be 
multiplied by a product of electron spin operators SiSj. Since an approximate azimuthal sym- 
metry around the magnetic (or atomic polarization) axis exists, one may integrate, or average, 
over the azimuthal angle <p of emitted photons, J Q 27r dep. This allows omitting various bilinears like 

oc cos (p sin if which all vanish after the azimutal angular integration. 
Thus, only two types of effective electron spin operators S±S T , S Z S Z , (S± = S x ± iS y ) contribute 
to rate calculation. 

For QED vertex we use two independent circular polarizations of h 7 = ±1, and write QED vertex 
amplitude d-E as circularly polarized electric field x ^ £ d + + ^ £ d_±-^d z with c = cos6*, s = sin 9 [9 

the photon emission angle from polarization axis) and d± = ^f(d x ±id y )/\/2. The squared amplitude 
didjEiE* is then bilinear function of the angle factor c, s. 

An example of parity violating (PV) term ending at the final state 1 3/2, —1/2) with emission of 
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This term gives the angular asymmetry oc sin 9 cos 9 for unpolarized (averaged over helicity /i 7 = ±1) 
photon emission, asymmetric with respect to 9 = tt/2. At the same time the same term gives circular 
polarization, when summed over the helicity weight Ylh ^7 X TS ^ e ^ leading to angular symmetric 
terms oc sin 9. 

Thus, unless other added terms cancel in total, two PV effects remain as oc cos r\ sin 77. 

Atomic wave functions we need for this Xe case are El dipole between the intermediate J = 1, 
and the final J=0 states, for which we use the measured decay rate j n f given by u 2 J2 po \ \d- E\ 2 /n = 
u;3 7™//^ra/ • The initial and the intermediate wave function overlap is due to the electron spin flip 
caused by the neutrino pair emission, hence having a large electronic orbital overlap. This part can 
be worked out approximately by using Clebsch-Gordon (CG) coefficients |10] alone. 

Spectrum rates and PV observables A general rate formula may be given under specified 
photon circular polarization h 1 = ±1, 



where c = cos#, s = sin#. Here ai,bi,ei,fi of mass dimension 2 are functions of variables 00 and 
elements of the neutrino mass matrix, containing product of 4 CG coefficients. Parity violating 
effects appear in and fi terms oc cos r\ sin r\. Parity conserving terms oc 6j,ej are multiplied by 



Without measurement of circular polarization h 11 there are components, symmetric one (even 
function of cos 9) relative to the atomic polarization axis, and asymmetric one (odd function of cos 9). 
The total rate oc 64/5 + 62/3 + 60 is an integral of symmetric functions. The angular asymmetry 
gives parity violating effect; for instance, the difference in the forward (F) and the backward (B) 
hemispheres at some particular angle 9 and n — 9oi0<9<7r/2. We plot in Figure 2 the angular 
asymmetry defined by the ratio of FB difference to FB sum, s^a^c 3 + a\c) / (64c 4 + 62c 2 + 60) , , taking 
several 9 of 7r/10, 7r/6, 7r/4, ir/3, 2tt/5 (with r\ = 7r/4). 

With measurement of circular polarization, the simplest P-odd quantity is the relative circular 
polarization as a function of angle summed at two angles 9 and tt — 9, given by s(f 2 c 2 + /o)/(64C 4 + 
62C 2 + 60) • We plot in Figure 3 their spectrum (taking 77 = vr/4) at the same angles as in Figure 2. 

Both of these PV quantities are of order unity, typically O[0.1], even far away from neutrino 
mass thresholds, since parity violation is intrinsic to the involved weak interaction not related to the 
finiteness of neutrino masses. The feature of large parity violation far from thresholds where the 
event rate is larger can be used for rejection of QED backgrounds and a positive identification of 
RNPE. The magnitude of PV quantities differ significantly at different emission angles. 

Parity-odd quantity thus computed is almost of equal magnitude in the Majorana and the Dirac 
cases, the difference being < 0[2x 10~ 4 ] in Xe case of a large energy difference ~ 8.4 eV. The Majorana 
vs Dirac distinction using P-even quantity is made easier by using atomic system of smaller energy 
difference close to neutrino pair masses such as 2777.3 (the largest pair mass); the rate difference may 
become several xlCT 2 for A < leV. 



d 2 T 




2 l c ^'| 2 (^ 4c4 + ^ 2 ° 2 + ^° + s ( a 3C 3 + aic) 



duod cos 9 



+/i 7 (e 3 c 3 + e lC + s{f 2 c 2 + / ) )) 
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cos 2 rj or sin 2 77. 
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Figure 2: Angular asymmetry of emitted photon at several angles 9 from oriented 131 Xe atoms in a 
metastable J = 2,F = 7/2 state as a function of photon energy u>. Assumed neutrino masses are 
lmeV, 10 meV, 50meV, with mixing angles, sin 6 12 = 0.3, sin 2 13 = 0.039 assuming the normal 
hierarchy [5J. 
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Figure 3: Circular polarization for the Xe metastable state at several angles from oriented atoms. 
Assumed parameters are the same as in Figure 2. 
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Macro-coherent RNPE Coherence of targets is important, both to enhance the event rate 
and to determine the threshold location. Ordinarily, the momentum conservation in elementary 
processes does not play important roles in atomic radiative decay, since atomic recoil is very small 
and difficult to measure. We however utilize the macro-coherence of target atoms, similar to, but 
more enhanced in rates than (hence called the macro-coherence), a single photon superradiance 

In macro- coherent phenomena more than two light particles such as 77 in two photon emission 
and ^UiUj in RNPE are involved in the final state of decay product. If the phases of emitted plane- 
wave particles of product amplitude e l ^ iki ' Xd from various target sites at Xd (distributed uniformly 
with a number density n) are coherently added, one has the macro-coherent rate in the large target 

volume (V) limit pj, T = nN j (lli^) (2tt) 4 5(A - £\ u;;) <5Q~\ A*) \M{h)\\ (N = nV), a result 
of coherence not restricted by the wavelength of emitted photon, quite unlike in the single photon 
superradiance [II] . Thus, the coherence prevails in a macroscopic target region, and one expects the 
momentum conservation, k + p± + P2 = for RNPE, and a large enhanced rate. 




(0{eV) 



Figure 4: RNPE spectrum rate for the 131 Xe metastable state in the unit of nNG 2 F , y n f / (2tt 2 A^f) of 
eq. (J3]) . The red-dotted threshold region including all six is expanded in the inset to exhibit threshold 
rises. Assumed parameters are the same as in Figure 2. 

Coherence dynamically develops, and the key for this is the seed of quantum fluctuation of fields 
or succesive spontaneous emission within the active target volume, as emphasized in superradiance 
phenomena [TTj . Moreover, an active external trigger by a pulsed laser may help much for the 
coherence development [13J. Analysis based on the Maxwell-Bloch equation [TTJ is expected to give 
strong angular correlation along the axis of exciting and triggering laser, an effect not evident in the 
present argument and ones in [3] as well. 

The total rate J dudT/cLo from the entire target roughly scales with nNG 2 F ^ n f / '(2ii 2 A 3 ^) times 
an integrated number. This is numerically ~4x 10 _9 (n/10 24 cm _3 )(A^/10 20 )(7 n //10 9 s~ 1 ) Hz times 
an integrated spectrum in the eV 4 unit. A 100 gr of solid 131 Xe gives the total rate of order ~ 3 x 10~ 5 
Hz, taking into account the enhancement due to the soliton formation crudely estimated in [5] to be 
O[10 4 ]. One needs a more detailed numerical integration of Maxwell-Boltzmann equation to give a 
more accurately determined rate. 
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The elementary rate from a single atom may be estimated by taking N = 1, and replacing the 
target number density n by the phase space density ~ A 3 of atomic decay, which is ~ wavelength -3 . 
Thus, the macro- coherent rate is enhanced by the factor of the ratio of target volume/wavelength 3 
compared to the elementary rate. In contrast the single photon superradiance has the enhancement 
factor of order, (linear dimension/wavelength) 2 [TT] . 

We show in Figure 4 the total spectrum rate, taking 77 = 7r/4. A global feature of the macro- 
coherent spectrum is the sharp rise from thresholds and a fast decrease towards small photon energy 
region, ultimately oc u 3 . The step around 4.1994 eV is due to (33) mass threshold, whose location, 
but not spectral shape, is insensitive to the angle factor 9 13 if sin 2 6*13 < 0.039 (the present limit). 
In order to locate thresholds one needs accuracy of photon energy resolution of order 10~ 3 or even 
less. This accuracy is not problematic, since one can trigger RNPE by laser whose energy resolution 
is far better. The fraction that threshold regions occupy is however not large; for instance the ratio 
of integrated rate of the inset part to the total in Figure 4 is ~ 5 x 10 -4 for all 6 threshold regions. 
We also worked out the case of inverted hierarchy of neutrino masses; the difference between the two 
cases is not overwhelming, but a sizable spectral shape difference exists between the two cases, the 
normal and the inverted hierarchies. In another word one can equally investigate and experimentally 
distinguish the normal and the inverted cases by this type of atomic experiments. 

In this Xe example it is necessary to obtain large high statistics data near threshold regions for 
the precision neutrino mass spectroscopy. This is due to that Xe level spacing ~ 8.4eV is somewhat 
large compared to the anticipated neutrino mass ~ 0.05eV in the assumed normal hierarchy case. 
Atoms of much smaller energy spacing are far better for the purpose of precision neutrino mass 
spectroscopy. 

To obtain a large enhanced macro-coherent rate, one needs a collection of atomic targets of 
0[Avogadro number] cm -3 , well isolated among themselves and weakly interacting with envirorn- 
ment. An example is a matrix implantation of a fraction of targets within solid matrix such as 
para-H 2 . 

In summary, we discussed two types of parity violating effects; angular asymmetry from oriented 
atoms and emergent circular polarization, as a means to positively verify weak RNPE process, which 
can be used to systematically determine the neutrino mass matrix. 
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